Spontaneous Isotropy Breaking: A Mechanism for CMB Multipole Alignments 



in 
o 
o 

(N 

Or 
<D 
00 

(N 



O 

m 

On 
O 

in 
o 

a: 

6 
to 

in 
c3 



Christopher Gordon, Wayne Hu, Dragan Huterer and Tom Crawford 
Kavli Institute for Cosmological Physics, Department of Astronomy and Astrophysics, 
and Enrico Fermi Institute, University of Chicago, Chicago IL 60637 

We introduce a class of models in which statistical isotropy is broken spontaneously in the CMB 
by a non-linear response to long-wavelength fluctuations in a mediating field. These fluctuations 
appear as a gradient locally and pick out a single preferred direction. The non-linear response 
imprints this direction in a range of multipole moments. We consider two manifestations of isotropy 
breaking: additive contributions and multiplicative modulation of the intrinsic anisotropy. Since 
WMAP exhibits an alignment of power deficits, an additive contribution is less likely to produce the 
observed alignments than the usual isotropic fluctuations, a fact which we illustrate with an explicit 
cosmological model of long-wavelength quintessence fluctuations. This problem applies to other 
models involving foregrounds or background anisotropy that seek to restore power to the CMB. 
Additive models that account directly for the observed power exacerbate the low power of the 
intrinsic fluctuations. Multiplicative models can overcome these difficulties. We construct a proof 
of principle model that significantly improves the likelihood and generates stronger alignments than 
WMAP in 30-45% of realizations. 



I. INTRODUCTION 

Data from the first year of WMAP has largely con- 
firmed the standard cosmological model based on nearly 
scale-invariant statistically homogeneous and isotropic 
fluctuations [1]. Yet interestingly the large angle deficit 
in temperature fluctuations versus the standard model, 
first seen by COBE [2], was both confirmed and sharp- 
ened. This result is unlikely at between 0.7% and 10%, 
depending on the assumptions, the choice of statistic, 
and the analysis [3-6]. Cosmological explanations for a 
deficit of power involve breaking scale-invariance initially, 
dynamically, or topologically. 

More recently, however, several other anomalies have 
been noted that point to a possible breaking of statistical 
isotropy on large angles. Statistically isotropy requires 
that all of the multipole moments of the CMB tempera- 
ture field be uncorrelated at the two point level. Yet the 
WMAP data exhibit curious correlations or alignments in 
the temperature field. In particular, the octopole of the 
CMB is highly planar and aligned with the quadrupole 
[7, 8]. Specifically the normalized angular momentum of 
the octopole is anomalously large at the 1/20 level and 
the angular momentum directions of the quadrupole and 
octopole are anomalously close at the 1/60 level. These 
and other alignments in the large-angle CMB have been 
further discussed and analyzed by Refs. [5, 9-13] and 
comprehensively studied in Ref. [14]. 

While the a posteriori nature of these statistics makes 
a strict interpretation of these probabilities difficult, the 
anomalies are clearly significant enough to merit serious 
attention. There are at least three classes of physical 
explanations: systematic, astrophysical, and cosmologi- 
cal. The WMAP instrument and time stream [15] have 
passed stringent tests that would make a systematic er- 
ror of this magnitude and type difficult to conceive. An 
astrophysical explanation involving foreground contami- 
nation is certainly possible. However known foregrounds 



are accounted for very successfully (e.g. [16, 17]) and the 
contamination would need to be of order the intrinsic 
CMB temperature anisotropy to account for the full ef- 
fect. Still foreground and systematic contamination may 
weaken the significance of at least some of the statistical 
anomalies typically at the expense of others [13, 14, 18]. 

The final possibility is that the anomalies have a 
cosmological explanation. Large-scale homogeneity and 
isotropy are clearly good assumptions for the universe as 
a whole. In this study, we examine the possibility that 
statistical isotropy is spontaneously broken in the fluc- 
tuations observed from a given spatial position and not 
in the fundamental theory. Correspondingly we will take 
the all-sky WMAP data [7] at face value and ignore the 
possibility of foreground and systematic contamination. 
We begin in §11 with the general mechanism and proceed 
in §111 and §IV to describe two different implementations 
of the mechanism. We discuss these results in §V. 



II. SPONTANEOUS ISOTROPY BREAKING 

We introduce a general mechanism for the breaking of 
isotropy in large angle fluctuations in §11 A. Under this 
mechanism statistical isotropy is preserved in the full the- 
ory but is spontaneously broken due to long-wavelength 
field fluctuations that appear as a gradient locally to the 
observer. A non-linear response to the field by CMB 
temperature fluctuations carries this preferred direction 
into a spectrum of multipoles. In §11 B we discuss the 
statistical treatment of the induced correlations between 
multipole moments used in the following sections. 



A. Non-linear Gradient Modulation 

There are two ingredients in our mechanism for break- 
ing the statistical isotropy. First, we require the pres- 
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ence of a field whose spatial fluctuations are dominated 
by long- wavelength contributions. Locally, around an ob- 
server, the structure in the field looks like a pure gradient. 
Given that a gradient must pick out a single direction, 
its presence provides the preferred axis for the breaking 
of isotropy. Because the long wavelength field is itself 
a random fluctuation, the underlying full theory retains 
statistical homogeneity and isotropy. Statistical isotropy 
is only apparently broken locally, and we refer to this as 
the spontaneous isotropy breaking. 

The second ingredient is a non-linear response of the 
CMB temperature perturbations to this field. If the re- 
sponse were purely linear, then the field gradient would 
produce a dipole in the CMB. If the non-linearity were 
small, the higher order anisotropy would be increasingly 
suppressed with multipolc number. A familiar example of 
this latter case is the Doppler effect due to the motion of 
the solar system through the CMB. The observed dipole 
is of course created by this motion. The effect on the 
quadrupole is a small correction to the true quadrupolc 
proportional to v 2 , while the corrections to the higher 
multipoles are entirely negligible [2]. Here we seek to 
break isotropy in at least the quadrupole and octopolc. 
Therefore the non-linearity must be substantial. 

Let us now parameterize a class of model that possesses 
these ingredients. A pure gradient projected across the 
sky along some axis x 3 has the structure 

G(n) oc (1 + Gin • x 3 ) = (1 + d cos 6) , (1) 

where n is the sky direction, Gi is the relative amplitude 
of the gradient and 9 is the polar angle from the preferred 
axis. A non-linear response to this gradient can be de- 
scribed as a general function of G(n), F[G(n)]. This 
function can in turn be represented by its spherical har- 
monic coefficients with the pole aligned in the direction 
of the gradient axis 

F[G(h)] = w(h)=Y, w c Y ^- ( 2 ) 

l 

The azimuthal symmetry around the preferred axis of a 
gradient limits the coefficients to m = 0. 

In general, the CMB may exhibit a non-linear response 
in two different ways. It may produce an additive ef- 
fect that is uncorrelated with the intrinsic anisotropy 
or it may produce a multiplicative modulation of this 
anisotropy. Furthermore, since the intrinsic anisotropy 
arises from a multitude of physical effects, it is possible 
that only one of them carries a non-linear response. 

Thus let us model the temperature fluctuation field of 
the CMB as 

T(n) = ^(n) + /[l + «;(n)]B(n), (3) 

where A(h) and B(h) are statistically isotropic Gaussian 
random fields. In the trivial case where B(h) = 1, the 
non-linear response is additive with an amplitude pro- 
portional to the parameter /. Alternately A and B may 
represent two different physical sources of the intrinsic 



anisotropy, e.g. the Sachs- Wolfe effect and the integrated 
Sachs- Wolfe effect. The non-linear modulation then pro- 
duces a multiplicative breaking of isotropy in one of the 
effects. 



B. Statistical Framework 

It is convenient to describe the statistical properties 
of the model in terms of the multipole moments of the 
fields. The observed temperature field of Eq. (3) becomes 

T(n) = 5^t/ m y/ m (n), ( 4 ) 

lm 

and similarly for the underlying isotropic fields 

lm 

B(n) = X^mWn). (5) 

lm 

The assumption of statistical isotropy for the underlying 
fields A and B requires that their covariance matrices 
satisfy 

(a em ae> m >) — 5w5 mrn >C( a , 

( a £mbl'm>) = SfjjSmm'Cg , 

(b*l m bl'm'} = SfJJ 5mm> Cj! b . (6) 

However statistical isotropy is not preserved in the ob- 
served temperature field T(n). 

Taking the multipolc moments of Eq. (3), the product 
w(n)_B(n) becomes a convolution 

ttm = aim + fblm + / E R^b^m (7) 

lA 

with a coupling matrix written in terms of Wigner 3j 
symbols 

M _ , 1W ,/ (21+ 1X24 + 1X212 + 1) 
R em - (-1) V ^ 

X { J { m -m) ^ ' ^ 
The ensemble average of the multipolc moments becomes 

(t'lmtl'm) = £>«' [Ci a + 2/G" h + f 2 Cg b ] 

+f E #&4*<%- ( 9 ) 

lll' t l2 

Here the ensemble average is over realizations of the fields 
A and B with a fixed field w. The field T(n) exhibits 
broken statistical isotropy with different variances in m 
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FIG. 1: The fiducial cosmological model and the WMAP power in 
i averaged over m ( "avg" ) and from azimuthally symmetric ( "m = 
0") components in the dipole frame. In addition to an overall deficit 
of power in the quadrupole, the WMAP quadrupole and octopole 
exhibit a further reduction of power along the dipole direction. The 
fiducial cosmological model is defined by f2 m = 0.27, h = 0.72, 
Q b h 2 = 0.024, <5 C = 5.07 X lO" 5 and r = 0.17. 



modes of a given I and a covariance between different I 
modes. 

For the simplest case of an additive effect B(n) = 1 
and is not a stochastic field. We can however consider 
its effective power spectrum to be C\ h = Att5(q. The 
ensemble average of Eq. (9) reduces to 



<5«<C £ aQ + f 2 w e w e >5 m0 . 



(10) 



where fwi is a deterministic contribution or mean to 
the Gaussian distribution of ttQ. A multiplicative model 
manifests more complex phenomenology. Here each wi 
correlates the observed multipoles across a range A£ — I. 
With a spectrum of we the maximum range of the corre- 
lation is determined by the highest non-negligible mode 
in the spectrum. We will explore the additive and mul- 
tiplicative versions in §111 and §IV respectively. 

A hallmark of this mechanism is that the modulation 
is azimuthally symmetric and hence couples harmonics 
with the same m in the preferred frame. Note however 
that in a general coordinate system, for example galactic 
coordinates, different m values will also be correlated. 

Specifying the change in the coordinate system with 
Euler angles by rotating around the z-axis by (j) € [0, 2tt], 
then around the new y-axis by angle 9 S [0, tt] and finally 
around the new z-axis by angle 7 € [0, 2tt], the harmonic 
coefficients transform as (see for example [19]) 



5 T I 1 



(11) 



where the Wigner rotation matrix is given by 

BLhh, 0, <j>) = e-"^d mrh (9)e~^ (12) 



TABLE I: WMAP data vs. fiducial model in the dipole vs maxi- 
mum angular momentum direction (see §IVB). N(/Cf d shows the 
relative contribution of noise assumed. 



I m Cf d ^K 2 ) N e /Cr \t lm \ 2 /C, 
dipole L„ 
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1 
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3 
1 
2 
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4 
1 
2 
3 
4 

5 
1 

2 
3 
4 
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1233 



577 



322 



0.005 



0.007 



0.009 



202 



0.011 



0.005 
0.013 
0.396 
0.179 
0.079 
0.426 
2.155 
3.929 
0.156 
0.052 
0.271 
0.406 
0.035 
0.376 
0.035 
6.827 
0.077 
0.364 



0.005 
0.003 
0.406 
0.315 
0.003 
0.029 
2.560 
1.641 
0.939 
0.508 
0.401 
0.180 
0.014 
0.811 
1.254 
2.777 
2.769 
0.078 



and 



^mm(^) — 
min{£-t-m,£+m} 



E 



fc-max{0,m+m} 



COS - 

2 



2k—m — r 



sm- 



2l+m+rh-2k 



x( 1)k+ i +m y/(e + m)\(l-m)\(l + my.(e-lRj\ 
k\(£ + m — k)l{£ + rh — k)l{k — m — fh)\ 

is a real function. Thus the ensemble average becomes 



(tirntt'm') - /Mlm^' 



)D"~D t -, 

' mm mm s 



(13) 



and exhibits completely broken isotropy. 

In Tab. I and Fig. 1, we compare the power in the ob- 
served multipole moments of WMAP in the dipole frame 
to Cf d , the power spectrum of a fiducial isotropic ACDM 
model. This model, which we use throughout the paper, 
has fl m = 1 - Q Q = 0.27, h = 0.72, Q b h 2 = 0.024, 
S c = 5.07 x 10~ 5 and r = 0.17. The WMAP data is 
taken to be the cleaned full-sky map of Tegmark et al. 
(hereafter TOH; [7]). Note that the m — components 
of the quadrupole and octopole are low both with respect 
to the average over m in the WMAP data and compared 
with the fiducial model. Thus the power deficit in the 
quadrupole and octopole are aligned. With broken sta- 
tistical isotropy it is in principle possible to explain these 
and other anomalies in the power distribution of the data. 

By employing the TOH cleaned map we are ignoring 
residual foregrounds and systematics. While these may 
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or may not weaken the significance of alignments in the 
data [5, 13, 14, 16, 20], they do not affect the ability of 
a cosmological model to produce them. The latter is the 
main topic we address in this work. 



III. ADDITIVE ISOTROPY BREAKING 

To illustrate the idea of spontaneous isotropy breaking 
we first take the simple case of an additive contribution to 
the anisotropy. Spontaneous isotropy breaking requires a 
non-linear response in the observed CMB temperature to 
an underlying gradient field. We begin with a motivating 
example of a known systematic effect, the non-linear re- 
sponse of a detector, in §111 A. This non-linear response 
modulates the CMB dipole into higher order anisotropy 
aligned with the dipole. We then apply this idea to a 
cosmological example involving dark energy density fluc- 
tuations in §111 B. 

The statistical anisotropy caused by additive isotropy 
breaking generally produces alignments of excess power. 
While alignments of power deficits, like those seen in the 
quadrupole and octopole of the WMAP data, are possi- 
ble we show in §111 C that they are statistically less likely 
than a chance alignment in a statistically isotropic field. 
This conclusion holds true for any mechanism which seeks 
to restore power and hence statistical isotropy in the 
quadrupole and octopole by removing a contaminant that 
is uncorrelated with the intrinsic anisotropy. 



A. Instrumental Example 

A concrete but non-cosmological means of additive 
isotropy breaking is provided by the non-linearity of a 
detector responding to temperature fluctuations. Sup- 
pose that the detector responds to a true temperature 
signal of G(n) as 



T dct (n) = /F[G(n)]=/£> 



G(n) 



/ 



(14) 



where ct\ — 1. Here / is an arbitrary normalization scale 
for the non-linearity of the response. If aii>i ^ then 
^det 7^ G and the observed temperature is a non-linear 
modulation of the true temperature. 

The dominant temperature signal in a differencing ex- 
periment is the dipole arising from our peculiar mo- 
tion, G(n) = Tdip cos 9 in the dipole frame, with Tdi p = 
3.35mK. The non-linear response has the general form of 
Eqn. (3). Taking / = T dip , 



Tdip 



ctiPi(cos#) + a 2 
\2 



^P 2 (cos0) + 1 (15) 



+a 3 



Ip 3 (cos6) + Ip 1 (cos0) 





+ ... 



where Pi = \jAir/(2l + l)Yg are the Legendre polyno- 
mials, and therefore 



U>2 



w 3 



4 Fk 

3V5 a2 + --' 



(16) 



where . . . represent contributions from i > 3. With a 2 — 
a 3 <~ 10~ 2 , the 10~ 3 dipole anisotropy is modulated into 
a 10~ 5 quadrupole and octopole anisotropy with m = 
in the dipole frame. 

Given the observed multipole moments t£ m , the intrin- 
sic anisotropy becomes 



aim = tlm - fwi5 m0 . 



(17) 



For multipole moments that are observed to be anoma- 
lously low, as is the case for the quadrupole and octopole 
moments of WMAP in the dipole direction (see Tab. I), 
ttm — * and |a^o | 2 - * f 2 wjS m0 . Thus the removal of 
the additive contribution restores power to the intrinsic 
sky and can bring it closer to a statistically isotropic dis- 
tribution. Note that this restoration involves a chance 
cancellation between a Gaussian random realization of 
ai m and a fixed additive contribution fwi. This will be 
important in addressing the probability of such an occur- 
rence in §111 C. 

Detector non-linearity is unlikely to be the actual 
source of the observed alignments. In this model, the 
nonlinearity of the response would be at the percent level 
for the 3.35mK dipole. If this model were extrapolated 
to brighter sources, it would imply that ~ IK sources 
completely saturate the response which is demonstra- 
bly not true for the WMAP instrument. Nonetheless, 
the non-linear detector response model provides a con- 
crete and conceptually simple illustration of the isotropy 
breaking mechanism. It has the added benefit that it 
naturally picks out the dipole direction as the preferred 
axis along which the temperature is modulated as seen 
in the WMAP data. 



B. Dark Energy Example 

Cosmological examples of isotropy breaking follow a 
similar form but are more complicated to calculate. Here 
one must combine the dynamics of how the modulated 
field produces temperature fluctuations in physical space 
along with the projection along the line of sight onto 
angular space. The benefit of a cosmological model is 
that the spectrum wi of the additive contributions have 
a well defined shape and hence the model in principle has 
predictive power, e.g. an additive correction at relatively 
high multipoles predicts a specific change to the low order 
multipoles and vice versa. 

As a concrete example, let us consider the Integrated 
Sachs- Wolfe (ISW) effect induced by spatial perturba- 
tions in the dark energy density from a quintessence field 



5 



Q. If the perturbations to the field Q are dominated by 
superhorizon scale fluctuations then within our horizon 
the field will look linear 



Qi =B + B lX3 



(18) 



where Bi are constants, so that the local gradient picks 
out a preferred axis x 3 . The non-linear modulation is 
provided by the scalar field potential. Taking it to be of 
the form 

V(Q) =Vq[1 + f cos(Q/M )] , (19) 

the field gradient manifests itself as 

V(Qi) = V a [l + f cos(k x 3 + 5)} , (20) 

where fco = Bi/Mq and S = Bq/Mq. This initial poten- 
tial energy fluctuation is then transferred onto the CMB 
through the ISW effect and the dynamics of the field. 

As the actual field gradients are from superhorizon 
wavelengths, the field will remain frozen provided its 
mass m is small compared to the Hubble parameter H: 



H 2 



d 2 V 



dQ 2 



1 

H 2 



: 3 



d 2 v 



dQ 2 



V 



<3/ 



M p i 

M 



where M p \ 
Thus, if 



(8nG) 1 / 2 is the reduced Planck mass. 



M pl ^ 1 ■ 



(21) 



the field will not roll and the density field then carries 
potential energy fluctuations that are frozen. Its Fourier 
components are 



PQ 



(k) 



p 5 (27r) 3 ,5(k-fc x 3 ) 

+ {-e- lS (2n) 3 S(k + k Q K 3 ). (22) 



For the moment, we consider the case where only the 
dark energy is initially perturbed. The comoving curva- 
ture perturbation, £, will evolve if there is a perturbation 
to the total pressure, pr- As only the dark energy is per- 
turbed, initially the Universe is effectively homogeneous, 
Ci = 0, and £ only begins to grow during dark energy 
domination when the only non-negligible components to 
the total energy density px are the dark energy density 
and the pressureless matter, p m . The dark energy is as- 
sumed to be sufficiently light that it is approximately 
frozen and so pq ~ —Pq = const. The evolution of the 
comoving curvature perturbation, in a flat Universe, is 
given by (see for example the Appendix of [21]) 



C = Ci 



a da' Sp T 



PT +PT 



a da' 8pQ 



(23) 




0.01 



FIG. 2: Dark energy based additive contributions wp. Shown are 
the first few multipoles of \wi/sg\ as a function of ho/Ho- Long 
wavelength perturbations exhibit a spectrum in Wg that is steeply 
falling in i and is modulated by parity considerations from S£ . 



Denoting the relative energy densities of dark energy and 
matter today as Oq = pq/ptU=i and Q m = p m /pr\a=i, 
the usual rcdshifting of pressureless matter can be writ- 
ten as 

_ Pq 

Pm — n „ 

a 6 Uq 

Substituting this into Eq. (23) gives 
a 3 Sp Q £l Q 



C 



3 Pq Q r , 



(24) 



The Newtonian gravitational potential can be expressed 
in terms of C as (see for example [21, 22]) 



T /, \ , H f a da' . 



5p T 



PT +PT 



, H f a da' , 



(25) 



Substituting Eqs. (23) and (22) into this gives 

*(k, a) » V(2ir) 3 tf(k-A:oi3)+V* (27r) 3 ,5(k+fc x 3 ) , (26) 
where 



da' ,3 
-a 



H(a') 



(27) 



When there are only dark energy perturbations, the 
temperature fluctuations in the direction n due to the 
potential ^ are given by the ISW effect. Relating these 
contributions to the form of an additive isotropy breaking 
term in Eq. (3) gives 



AT(n) 



fw(h) = J 2*'(x = Dh,a)d\na, 



(28) 
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where D is the comoving distance along the line of sight 
and primes are derivatives with respect to In a. The mul- 
tipole moments are given by 

r at 

fw t = j dnY/ m (n) — (n). (29) 

With the Rayleigh expansion of a plane wave 

exp(ik • x) = Wj,(fcL>)r; m (k)y, m (n) , (30) 
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they can be written as 

d 3 k 



din a 



(2tt) 3 



W J ,(i : fl)y; m (k)2$'(k,«). 



(31) 

Taking the form of the Newtonian gravitational potential 
above, we get 

w e = --2- St y/4Tr(2e+l) / d In a Je (k D)I(a)a 3 , 

s t = cos^(-l) £ / 2 (5| + sin5(-l)( £+1 )/ 2 ^ , (32) 

where Sg — 1 if I is even and if £ is odd, and vice versa 
for 5?. Here we have defined 



1(a) 



3a 3 d In < 



H(a) 



da 1 



/3 



(33) 



Given that the CMB response to a dark energy density 
perturbation is still linear, restoring adiabatic fluctua- 
tions is a simple matter of adding in the fiducial model 
contributions as an uncorrelated contribution to the tem- 
perature field. 

In Fig. 2, we show the amplitude of the first few mul- 
tipoles of we/se as a function of ko/Ho in the fiducial 
cosmological model. There are two general features in 
we that go beyond the specifics of this model. First, a 
spatial modulation projects onto an angular modulation 
with a weight given by the spherical Bessel function je- 
For a supcrhorizon fluctuation 



kn kn 

k D = -r^HoD ~ JL < 1 , 

f no 



(34) 



and so je oc (fco/i/o) • Thus a superhorizon scale mod- 
ulation typically predicts a sharply falling spectrum in £ 
for the additive contribution. The modulation scale ko 
should not be confused with the field gradient scale which 
is by assumption always superhorizon. 

The symmetry factor se provides an interesting modu- 
lation of this result, and this is the second general feature 
of the model. The phase of the plane wave 8 determines 
its symmetry under a parity transformation n — ► — n 
or X3 — > X3. Because spherical harmonics have definite 
parity (— l) e , the plane wave modulation affects even and 
odd multipole moments differently. In particular choos- 
ing a phase S — > 0, eliminates the contributions from all 
odd multipoles. 
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FIG. 3: Ax 2 relative to the fiducial model for the additive model 
parameter space of phase and wavelength (<5, fco / Ho ) with / chosen 
to minimize the Ax 2 . For the multipoles I = 4 — 10, the maximum 
improvement is Ax|_ 10 = —5.4 (shaded contours). These models 
restore power to the quadrupole and octopole and thus £ = 2,3 add 
a disfavored positive contribution through Ax|_ 3 (curves). 
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FIG. 4: Restoration of power in the dark energy additive model. 
Taking a phase of <5 = — 0.173-7T that intersects the minimum 
Ax|_ 10 from Fig. 3, models with wavelengths ko/Ho around the 
minimum restore power to the quadrupole and octopole (dashed 
line vs points) and bring them closer to the fiducial model (curve). 



C. Power Excess and Restoration 

Given the dark energy based model of the previous 
section, we can test the WMAP data for additive contri- 
butions that break statistical isotropy. We set the pre- 
ferred axis for isotropy breaking to be in the direction of 
the dipole based on the observed pattern of low multi- 
pole alignments (e.g. [8, 10, 11]). In the dipole frame, 
the observed power in te m is given as a fraction of the 
power in the fiducial model of Fig. 1 in Tab. I. The un- 
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even distribution of power across m for the quadrupolc 
and octopole quantifies the alignment problem. 

To determine whether the WMAP data favor an ad- 
ditive breaking of statistical isotropy, we maximize the 
likelihood with respect to model parameters. Since the 
additive contribution is deterministic and the underly- 
ing temperature field is assumed to be drawn from the 
fiducial model this amounts to minimizing the x 2 

^max 

X 2 = Y1 ( Ce + A ^)~ 1 (*lm - fw e 8 m0 )(t em - fw e 8 m0 ) . 

. . (35) 

Here we have included a crude estimate of the noise in the 
WMAP data from their quoted power spectrum errors in 
the absence of sample variance. These are listed in Tab. I. 
Note that our results here are insensitive to the specific 
representation of the noise as it is much smaller than the 
rms signal of the fiducial model. 

The benefit of having a cosmological model for the ad- 
ditive contributions is that the functional form of the 
contributions is determined through wg. Evidence for an 
alignment of excess power in one range of multipoles pre- 
dicts the correction that must be applied to another. We 
therefore minimize the \ 2 for 4 < I < 10 with respect 
to the parameters (/, kg/Ho^) and then examine what 
it implies for the quadrupole and octopole. The best-fit 
solution is k /H = 4.8,6= -0.173tt,/ = 1.16 • 1(T 4 and 
has improvement of A% 2 = —5.4 relative to the fiducial 
model with no additive contribution. In fact, there is a 
wide range of models with a similar level of improvement 
(A% 2 <~ —5) since the two parameters k /H n and S affect 
the I dependence of the contributions similarly. We show 
in Fig. 3 the A% 2 contours for I = 4 - 10 in the (k /H , 
S) space for the best-fit amplitude /. About two thirds 
of the improvement in \ 2 comes from £ = 4 which has a 
~ 2er excess in power compared with the fiducial model 
(see Fig. 1 and Tab. I). 

Given this range of models that improve the A% 2 for 
I = 4 — 10, it is interesting to see what they imply for the 
quadrupole and octopole. In Fig. 4, we show the intrin- 
sic WMAP power spectrum (i.e. one with the additive 
contribution removed) 

2i + i \ tlm ~ / w ^™o| 2 (36) 

m 

for 5 = — 0.1737T and a range of ko/Hg that includes the 
best fit model. Models around the best fit restore power 
to the quadrupole and octopole to bring them closer to 
the fiducial model. 

There are however a number of problems that make 
this explanation of the alignment problem unsatisfactory. 
First, although the removal of the excess power in I = 4, 
m = brings the distribution in m of the power closer to 
flat, the total power is well under the fiducial model. Our 
additive model is limited in that it can only affect m = 
in the preferred frame. Similarly the quadrupole and 
octopole multipole moments are low not only in m = 



but also m = ±1 (see Tab. I) which cannot be altered in 
our additive model. 

The most important flaw however is that to restore 
power to the quadrupole and octopole one must assume 
that the additive contribution nearly exactly cancels the 
intrinsic temperature fluctuation. To assess whether that 
is a likely occurrence, consider the A% 2 contributed by 
1 = 2 and 3. Around the minimum for t = 4 — 10, it 
receives a positive contribution of A%2-3 between 1 and 
4 relative to the fiducial model (see Fig. 3). A chance 
cancellation of power is less likely than a realization of 
anomalously low power. For example, in random real- 
izations of the fiducial model without the additive cor- 
rection we find 2.3% of the realizations have a lower av- 
erage power in the quadrupole moments than WMAP. 
With the best fit additive correction this fraction drops 
to 0.8%. 

Note that quadrupole and octopole problem remains 
even had the improvement from Ax 2 _ 10 been signifi- 
cantly larger. It also remains if the figure of merit were 
the alignment of the quadrupole and octopole and not 
X 2 . The reason is that the Gaussian distribution for the 
intrinsic sky in a^ m peaks at zero. The probability that 
a realization of falls within 5aio ~ \teo\ of the addi- 
tive contribution W( is uniformly lower than the prob- 
ability that it fell by chance at a \aeo\ < 8a^ without 
wg. The same argument holds for alignments of power 
deficits since they are formed from the joint probability 
of low power across a range of I. 

These conclusions hold for any additive model where 
a template contamination by chance cancels the intrin- 
sic fluctuations. This includes explanations of the align- 
ments involving foregrounds and broken isotropy in the 
background metric as in Bianchi models [23]. A can- 
cellation model can at best be said not to decrease the 
probability of alignments of low power multipoles by a 
substantial amount. It can never increase the probabil- 
ity. 

There is one remaining possibility for additive models. 
An additive contribution that is not azimuthally sym- 
metric around the preferred axis might alone account for 
essentially all of the observed quadrupole and octopole. 
For example, foreground contamination could in principle 
supply the observed power. Unfortunately, such an ex- 
planation of alignments would then exacerbate the large 
angle power deficit problem. 



IV. MULTIPLICATIVE ISOTROPY BREAKING 

A multiplicative model for breaking statistical isotropy 
eliminates the two fundamental flaws of the additive 
model. Firstly, though the modulation is still az- 
imuthally symmetric it affects multipole moments of all 
m values. Secondly, the change in the observed multipole 
moments is correlated with the intrinsic anisotropy and 
hence can explain alignments of power deficits. Together 
these features can explain the planarity and alignment of 
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FIG. 5: Multiplicative model with a quadrupolar modulation. 
Shown is the fractional change or transfer of power from an in- 
trinsic multipole (£, m) to an observed (£',m') or p% m j_^/ m ,y 
For a modulation amplitude —5 < W2 < —10, power in i = 2, 3 and 
m = 0, 1 is suppressed while that in (£,m) = (2,2), (3,3), (4,0) 
and (5, 3) is enhanced. These features are seen in the WMAP 
TOH data of Tab. I and enhance the planarity and alignment of 
the quadrupole and octopolc. 

the quadrupole and octopole seen in WMAP. 

Based on the instrumental example of a calibra- 
tion error that varies across the sky (§IVA), we con- 
struct a proof of principle model (§IVB) that solves 
the quadrupole and octopole mutual alignment problem 
(§IVC). 

A. Instrumental Example 

An instrument whose calibration varies across the sky 
will modulate the intrinsic anisotropy -B(n) as 

T(n) = f[l + w(n)]B(n), (37) 

where / is the sky averaged calibration and w is the fluc- 
tuation in the calibration. A calibration fluctuation can 
transfer power from one multipole moment to another as 
well as alter the power in a given multipole. 

Consider for a calibration fluctuation that is purely 
quadrupolar in nature 

w(h) =w 2 Y 20 {n). (38) 



This form is motivated by the need to suppress fluctu- 
ations along both poles of the preferred axis. It also 
relates multipoles of A£ = 2 by Eq. (9) (see also [24] for 
a dipolar example) which as we will see has interesting 
implications for alignments and parity anomalies beyond 
the quadrupole and octopole [11, 25]. 

The effect of this calibration modulation is an £ and 
m dependent alteration of the amplitude of the temper- 
ature field as well as a transfer of fluctuations between 
A£ — 2 modes of the same m. From Eq. (7), the contri- 
bution from an intrinsic fluctuation bt> m i to an observed 
multipole moment tg m is given by 

fb r rn, (39) 

where 

P(t> ,m')->(l,m) = ($££'+ Rhn^J dmm> ■ (40) 

These square of several of these p's are plotted as a func- 
tion of u>2 in Fig. 5. For certain values of 1112, the modu- 
lation can suppress specific multipoles of the quadrupole 
and octopole to zero. Interestingly for — 10 < W2 < —5, 
the suppression is efficient for m = 0, ±1. Moreover for 
this range, the power in I = 2, m = ±2 and £ = 3, 
m = ±3 is actually enhanced. These properties are ex- 
actly what is needed to explain the observed power dis- 
tribution of the WMAP sky. 

Furthermore, a modulation in this regime also en- 
hances power in £ = 4, m — from intrinsic power in 
I = 2, m = and I = 5, m = 3 from i = 3, m = 3. The 
WMAP multipoles (4, 0) and (5, 3) have excess power in 
the dipole direction. Finally the modulation correlates 
the observed signal between these pairs of multipoles. 
The correlation for £ — 2, m = is not very relevant since 
its observed value is consistent with noise (see Tab. I). 
For £ — 3, to = 3 it is of the right sign to account for 
enhanced fluctuations in both multipoles as seen in the 
WMAP data [11]. 

Unfortunately a simple calibration error cannot lit- 
erally be the cause of alignments. In the discussion 
above we have only considered the modulation of intrinsic 
power in the quadrupole and octopole. A calibration er- 
ror would modulate all of the intrinsic anisotropy. In this 
case terms such as P(4,o)-> (2.0) would generate observed 
quadrupole power out of intrinsic £ = 4 fluctuations and 
spoil the ability to lower its power. 



B. Proof of Principle Example 

At least in principle, a cosmological version of the 
calibration mechanism has the ability to modulate only 
the quadrupole and octopole. Large scale temperature 
anisotropies can arise from both the Sachs- Wolfe and in- 
tegrated Sachs- Wolfe effect and each can carry extra de- 
grees of freedom from isocurvature fluctuations or multi- 
ple sources of curvature fluctuations. It is possible then 
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FIG. 6: ^Xcff = — 2Aln£ for the multiplicative model with the 
dipole as the preferred axis. Models with W2 ~ —5 and M3 £S 0.8 
have Ax^g ~ — 16 and best match the dipole frame data given that 
the most anomalous power deficit is in £ = 2, m = (see Fig. 5). 
The amplitude / is set to minimize A^g for each model. A range 
of models around the maximum improve the Ax?<j and most of the 
improvement is due to the breaking of statistical isotropy through 
the single parameter W2- The dashed box highlights models which 
best produce angular momentum alignments in Fig. 7. 



that the presence of a long wavelength field gradient mod- 
ulates only one of these effects. 

Constructing a fully physical model is beyond the 
scope of this paper. Here we seek only a proof of prin- 
ciple that a multiplicative modulation of the intrinsic 
anisotropy can explain the WMAP anomalies. This proof 
of principle should help to guide the construction of an 
physical model that satisfies the observational require- 
ments. 

Let us take the general form for the modulation Eq. (3) 
where the anisotropy is constructed out of two contribu- 
tions A(h) and B(ri), only one of which is modulated. 
Based on the discussion of calibration error in the pre- 
vious section, we parameterize their underlying power 
spectra as 



Cf a = (1 
Cf = 0, 



M e )C, 



fid 

e ' 



M e C, 



lid 



(41) 



Here Cf d is the fiducial cosmological model specified in 
§111 B. We will further take the modulation to be purely 
quadrupolar in nature as given by Eq. (38). As discussed 
in §111 B, a nearly pure quadrupolar modulation can be 
achieved in physical space by the combination of a long 
wavelength spatial modulation whose phase corresponds 
to a nearly even parity when projected on the sky. 

The multiplicative modulation model is then specified 
by three quantities: w-i which controls the relative am- 
plitude of the modulation in the field B(h), f which de- 
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FIG. 7: Percentage of multiplicative models that exhibit a higher 
aligned angular momentum statistic L 2 for the quadrupole and 
octopole than WMAP. Here / is chosen to minimize the Axj^j of 
the data in the dipole frame. The percentage peaks at 45% for 
W2 ~ —7 and M3 = 1 and is to be compared with 0.22% for the 
fiducial model. 



scribes the relative amplitude of A(n) and B(n) contribu- 
tions and Mi which describes their relative contribution 
to each multipole moment. Given the discussion above, 
we seek a model where the intrinsic anisotropy is maxi- 
mally modulated for the quadrupole and minimally mod- 
ulated beyond the octopole. Note that the model is con- 
structed so that in the absence of modulation (w2 = 0) 
Ce = Cf d for / = 1. Hence for simplicity we set M% = 1 
and Mi>i = 0. The proof of principle model is then 
defined by three parameters u>2, M3, and /. 



C. Power Deficits and Alignments 

We now assess whether the WMAP data favor such a 
multiplicative modulation. We again begin by consider- 
ing the likelihood 



£(t|C) 



1 



(2 7 r) iV / 2 VdetC 
exp[-Xeff 



exp 



1 



t+c- 



-V2] 



(42) 



where t is an N element vector of the WMAP TOH mul- 
tipole moments ti rn in the dipole direction. C is the sum 
of the signal covariance matrix in Eq. (9) and a diagonal 
noise covariance given by Ne in Tab. I. The signal co- 
variance matrix is parameterized by W2, M3 and / [see 
Eqs. (9) and (41)]. It is important here to include Ne 
since some of the multipole moments are consistent with 
noise and the model has the ability to lower the predicted 
power in certain multipoles to zero. On the other hand 
an accurate quantification of noise is not crucial. 
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In Fig. 6, we show the improvement in A%^ ff relative 
to the fiducial model as a function of (w2,M 3 ) for the 
value of / that maximizes the likelihood. The maxi- 
mum likelihood value at W2 = —5.0 and M3 = 0.82 
has a AXctf — —16- Even more striking, most of the 
gain comes through w% . The best unmodulated or sta- 
tistically isotropic model (restricted to W2 = 0) has only 
an improvement of &XcS ~ —5. Note that in this case 
maximization over (/, M3) corresponds to the a maxi- 
mization over (C^jCa). The maximum at W2 = —5.0 is 
driven by the most anomalously low multipole moment 
in the dipole direction I — 2, m — as shown in Fig. 5. 
Note however that a wide range of values surrounding 
the maximum would also improve the fit. 

To test the quadrupole-octopole alignment, we take the 
normalized angular momentum [8] as generalized by Ref. 
[14] 



intrinsic 



The statistic 



p YL=-e\ a z™\ 2 



(43) 



(44) 



maximized over direction of the preferred axis captures 
both the alignment of the quadrupole and octopole and 
the planar nature of the octopole. Using Eqs. (13) and 
(43) and maximizing over direction, we find the best fit 
axis to be 



I = 247?3, b = 59?7 



(45) 



in galactic coordinates. This is close to but differs from 
the direction of the dipole (I = 263?85 ±0°l,b = 48?25 ± 
0?04) [1]. For WMAP L^ MAP = 0.96. 

To assess whether the multiplicative modulation im- 
proves the probability of alignments, we draw Monte 
Carlo realization of the models. With the fiducial model 
only 0.22% of 10 4 realizations had a L 2 > ^wmap m an Y 
direction. For the multiplicative modulation models, the 
fraction of 10 3 realizations is shown in Fig. 7. Here / is 
fixed to minimize Ax^ ff for each {w2,M%) in the dipole 
direction 

This fraction is maximized at u>2 = — 7 and A/3 = 1 
at 45%, an increase of a factor of ~ 200 over the fiducial 
model. From Fig. 5, we see that W2 = — 7 achieves the 
best balance of reducing the power in all moments of the 
quadrupole and octopole except m = ±£. Hence such a 
model would be expected to produce the highest normal- 
ized angular momentum. In Fig. 8 we show an example 
of one of these realizations (with L 2 = 0.99) that illus- 
trates how the modulation takes an intrinsically isotropic 
sky into one that exhibits the alignments of WMAP. 

Moreover, we have verified that models that solve the 
alignment problems do so irrespective of which test is 
used to quantify them. For example, the aforementioned 
model with u>2 — — 7 and M3 = 1 makes the mutual 
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FIG. 8: A realization of the multiplicative model where the 
quadrupole (left column) and octopole (right column) exhibit an 
alignment similar to WMAP. First row: intrinsic (unmodulated) 
sky from a Gaussian random isotropic realization. Second row (sin- 
gle column): the quadrupolar modulation oc — [1 — 7l2o(ri)] in the 
dipole direction. Third row: the modulated sky of the observed 
CMB (with M z = 1). Fourth row: WMAP full-sky quadrupole 
and octopole. 



closeness of the quadrupole and octopole area vectors 
[10, 14] unlikely at 37% level, compared to 0.12% level 
without modulation. 

The alignment statistic also favors high values of M3 
since M2 = M3 = 1 implies that all contributions to the 
quadrupole and octopole are modulated. The x^s statis- 
tic slightly disfavors very high values since in the dipole 
direction there remains a substantial amount of power in 
£ = 3, m = ±1,±2 which is better fit with an uncorre- 
lated isotropic component (see Tab. I). Nonetheless there 
is a range of models around W2 = —7 and M3 = 0.95 
that both improve xis anc ^ produce a fraction of more 
extremely aligned skies of > 10%. 

In fact, even better models can be found if so desired. 
Choosing the preferred axis to be exactly the dipole for 
parameter fitting is not necessary. In fact in the i^ax di" 
rection of WMAP (see Eq. (45)) the planarity of the ob- 
served quadrupole and octopole is maximized and hence 
has particularly low m = ±l, m = ±2 components (see 
Tab. I). Fitting parameters of the model to the data 
with this as the preferred axis allows a wider range of 
solutions (see Fig. 9). Now models with more negative 
W2 are favored to reduce the power in these components 
and higher values of M3 are allowed. With the maximum 
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FIG. 9: Same as Fig. 6 except replacing the dipole as the pre- 
ferred axis with that of L^ ax . By maximizing the planarity of the 
quadrupole and octopole |m| < I modes are minimized allowing for 
smaller u>2 and higher M3. 
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FIG. 10: Same as in Fig. 7 except with / given by the mini- 
mum Ax^g in the L^ax direction. Alignment requirements remain 
largely unchanged from the dipole direction but maximum percent- 
ages are more compatible with substantial improvements in AxJ ff 
from Fig. 9. 



likelihood / values with this preferred axis, W2 = —7 and 
M 3 w 1 still produces more extreme alignments in ~ 45% 
of the Monte Carlos. Moreover models exist with XcS im ~ 
provements of more than —10 which produce alignments 
more extreme than WMAP in > 30% of the Monte Carlos 
(see Fig. 10). 

The change in Ax^ ff between the dipole and i^ ax di- 
rections and more generally its sensitivity to direction 
indicates that one should be careful in interpreting this 
improvement probabilistically. A A%^ ff w — 16 likelihood 



improvement naively suggests a "4-cr" result. However 
since the choice of the i^ ax direction maximizes the pla- 
narity of the quadrupole and octopole, there will be some 
preference for modulation even for realizations of a sta- 
tistically isotropic sky. 

To address this issue, we fit the model parameters in 
the 10 4 realizations of the isotropic fiducial model in the 
L^ ax direction of each run. We find that the WMAP 
measured sky shows a greater improvement in the like- 
lihood than 98.8% of the Gaussian random, isotropic 
skies indicating that the multiplicative model is indeed 
strongly preferred by the WMAP observations. 

Since this exercise is simply a proof of principle with 
somewhat ad hoc parameters, we do not pursue finding 
the best model parameters further. The important con- 
clusion is that models exist where a multiplicative modu- 
lation of large angle anisotropies in a quadrupolar pattern 
can make the WMAP alignments a likely occurrence. 



V. DISCUSSION 

We have introduced a general class of models that 
break statistical isotropy in CMB temperature fluctua- 
tions spontaneously. The key elements are a non-linear 
response in the CMB to long wavelength fluctuations in 
a mediating field. Long wavelength fluctuations pick a 
direction through their local gradient and a non-linear 
response allows this direction to be imprinted in a range 
of multipole moments. 

We have examined two subclasses of models based on 
whether isotropy is broken through an additive contri- 
bution to the anisotropy or a multiplicative one. The 
two classes have analogues in the instrumental domain 
and are related to saturation in detectors and an angular 
drift in calibration. 

As a cosmological example of the additive class, we 
considered gradients in a quintessence field that are non- 
linearly mapped through that field's potential into dark 
energy density perturbations. The CMB then responds 
linearly to these density perturbations through the inte- 
grated Sachs- Wolfe effect. The existence of such contri- 
butions allows for the underlying CMB to be isotropic. 
Nonetheless, additive contributions cannot be said to 
solve the alignment problem in WMAP. Since the anoma- 
lies in the quadrupole and octopole are related to a deficit 
of power near the dipole direction, an additive correction 
can only explain the observations if the intrinsic fluctua- 
tions cancel them in both. This is statistically less likely 
than having an aligned deficit of power by chance without 
isotropy breaking. 

This fundamental problem exists for all explanations of 
the alignments which seek to restore power to the CMB 
including those involving foregrounds or other uncorre- 
lated templates. Additionally, an additive contribution 
based on an azimuthally symmetric modulation of a gra- 
dient field can only affect m = multipoles around the 
preferred axis. The only remaining possibility is that 
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the additive contributions themselves account for most 
of the power in the quadrupolc and octopolc. All of the 
moments of the intrinsic quadrupole and octopole would 
then be even lower than those observed, exacerbating the 
power deficit problem. 

These difficulties can be overcome by a model which 
breaks isotropy multiplicatively. Here we have demon- 
strated a proof of principle that the alignments and power 
distribution of WMAP can be made likely with such 
a mechanism. By modulating the intrinsic quadrupolc 
and octopole in a quadrupolar fashion, the likelihood of 
the WMAP data can be increased by a factor of e 16 / 2 
and separately, the probability of obtaining a sky with 
a higher angular momentum statistic is increased by 200 
to <~ 45%. The angular momentum statistic is a quan- 
tification of both the planarity of the octopole and its 
alignment with the quadrupole. Note that the likelihood 
increase must be considered in the context of the direc- 
tional maximization of the angular momentum. The fidu- 
cial isotropic model would falsely produce greater than 
e i6/2 j m p rovemen t s j n ^ i% f realizations. 

In addition, the model predicts correlated power ex- 
cess power in £ = 3, m = 3 and I = 5, m = 3 and more 
generally the A£ = 2 parity pattern of anomalies in the 
data [25]. There is a range of parameters that jointly in- 
crease the likelihood by e 10 / 2 and the angular momentum 
probability to ~ 30%. 

There are also classes of anomalies that are not ex- 
plained by this version of spontaneous isotropy breaking. 
Solar system related alignments [10, 14] and the asym- 
metry in power between the northern and southern eclip- 
tic ecliptic hemisphere [26] both fall into this category. 
While isotropy breaking involving a dipolar modulation 
[24] along the ecliptic could potentially explain the latter 
feature, it cannot simultaneously explain anomalies along 
the dipole direction. Furthermore, there is no fundamen- 
tal link between the dipole of our peculiar motion and the 



preferred axis for symmetry breaking. That coincidence 
remains in the multiplicative model. 

We emphasize that this proof of principle is not a phys- 
ical model and challenges will have to be overcome in con- 
structing one. The level of multiplicative modulation pre- 
ferred is fairly extreme and cannot be considered a per- 
turbation. It must affect essentially all of the quadrupolc 
and octopole in spite of the fact that in the fiducial cos- 
mology there are two physically distinct sources to these 
fluctuations: the Sachs- Wolfe and ISW effects. Thus a 
physical solution will have to modify the power spectrum 
of these contributions as well. On the other hand the 
modulation must not continue to intrinsic fluctuations 
in multipoles much higher than the octopole. Thus two 
distinct sources of anisotropy are required. 

Constructing a physical mechanism that satisfies all of 
these requirements will clearly be challenging and is be- 
yond the scope of this paper. Likewise the significance of 
the alignments and hence the model improvements may 
be degraded by foreground and systematic contamina- 
tion not accounted for in the all-sky cleaned maps, at 
the expense of exacerbating the power deficit problems. 
Nevertheless, this proof of principle model can serve as a 
guide both to building a physical model and to searching 
for alternate explanations of the anomalies. 
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